Chapter 16 



Transporting densities 



Paulina: I'll draw the curtain: 

My lord's almost so far transported that 

He'll think anon it lives. 

— W. Shakespeare: The Winter's Tale 

(P. Cvitanovic, R. Artuso, L. Rondoni, and E.A. Spiegel) 

In chapters 2, 3, 7 and 8 we learned how to track an individual trajectory, and 
saw that such a trajectory can be very complicated. In chapter 4 we stud- 
ied a small neighborhood of a trajectory and learned that such neighborhood 
can grow exponentially with time, making the concept of tracking an individual 
trajectory for long times a purely mathematical idealization. 

While the trajectory of an individual representative point may be highly con- 
voluted, as we shall see, the density of these points might evolve in a manner that 
is relatively smooth. The evolution of the density of representative points is for 
this reason (and other that will emerge in due course) of great interest. So are 
the behaviors of other properties carried by the evolving swarm of representative 
points. 

We shall now show that the global evolution of the density of representative 
points is conveniently formulated in terms of linear action of evolution operators. 
We shall also show that the important, long-time "natural" invariant densities are 
unspeakably unfriendly and essentially uncomputable everywhere singular func- 
tions with support on fractal sets. Hence, in chapter 17 we rethink what is it that 
the theory needs to predict ("expectation values" of "observables"), relate these 
to the eigenvalues of evolution operators, and in chapters 18 to 20 show how to 
compute these without ever having to compute a "natural" invariant density po. 
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Figure 16.1: (a) First level of partitioning: A 
coarse partition of M into regions Mo, M\, and 
Mi- (b) n = 1 level of partitioning: A refinement 
of the above partition, with each region A4, subdi- 
vided into Mm, Ma, and Ma- 



(a) 




(b) 



16.1 Measures 



Do I then measure, O my God, and know not what I mea- 
sure? 

— St. Augustine, The confessions of Saint Augustine 



A fundamental concept in the description of dynamics of a chaotic system is that 
of measure, which we denote by d/u(x) = p{x)dx. An intuitive way to define and 
construct a physically meaningful measure is by a process of coarse-graining. 
Consider a sequence 1, 2, n, ... of increasingly refined partitions of state space, 
figure 16.1, into regions At, defined by the characteristic function 



A coarse-grained measure is obtained by assigning the "mass," or the fraction of 
trajectories contained in the ith region A1, c At at the nth level of partitioning of 
the state space: 



The function p(x) = p(x, t) denotes the density of representative points in state 
space at time t. This density can be (and in chaotic dynamics, often is) an ar- 
bitrarily ugly function, and it may display remarkable singularities; for instance, 
there may exist directions along which the measure is singular with respect to the 
Lebesgue measure (namely the uniform measure on the state space). We shall 
assume that the measure is normalized 



where the sum is over subregions i at the nth level of partitioning. The infinites- 
imal measure p(x) dx can be thought of as an infinitely refined partition limit of 




(16.1) 




(16.2) 




(16.3) 
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Figure 16.2: The evolution rule /'can be used to map 
a region M; of the state space into the region f'(Mj). 




A//; = |At/| p(jt/) , where |Al/| is the volume of subregion At; and x\ e At/; also 
p{x) is normalized 



f Jxp(x) - 1 . (16.4) 

Here |At/| is the volume of region At;, and all |Al/| — > as n — > oo. 

So far, any arbitrary sequence of partitions will do. What are intelligent ways 
of partitioning state space? We already know the answer from chapter 11, but let 
us anyway have another look at this, in order to develop some intuition about how 
the dynamics transports densities. chapter 1 1 



16.2 Perron-Frobenius operator 



Given a density, the question arises as to what it might evolve into with time. 
Consider a swarm of representative points making up the measure contained in a 
region At/ at time t = 0. As the flow evolves, this region is carried into /'(At/), 
as in figure 16.2. No trajectory is created or destroyed, so the conservation of 
representative points requires that 



I dxp(x,t)= J dx p(xo,0). 
Jf'(Mi) JMi 



Transform the integration variable in the expression on the left hand side to the 
initial points xq = f~\x), 

I dx p(f t (xo),t)\detJ , (xo)\= I dx p(x ,0). 

The density changes with time as the inverse of the Jacobian (4.28) 

p(x,t)= x = f\x Q ), (16.5) 

|det J'(x Q )\ 

which makes sense: the density varies inversely with the infinitesimal volume 
occupied by the trajectories of the flow. 
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The relation (16.5) is linear in p, so the manner in which a flow transports 
densities may be recast into the language of operators, by writing exercise 16.1 

p{x,t) = (£! op){x)= f dx s(x-f(x ))p(x ,0). (16.6) 
Jm 



Let us check this formula. As long as the zero is not smack on the border of dM, 
integrating Dirac delta functions is easy: J M dx 5{x) = 1 if € At, zero otherwise. 
The integral over a 1 -dimensional Dirac delta function picks up the Jacobian of its 
argument evaluated at all of its zeros: 



(x—x )h'(x ) 



f 



dx6(h(x)) 




y — 



{x:h(x)=0\ 



(16.7) 



and in d dimensions the denominator is replaced by 



dx S(h(x)) 



-u 



Mj 



dxS(h(x)) = 



det 



dh(Xj) 
dx 



(16.8) 



where Mj is any open neighborhood that contains the single Xj zero of h. Now 

you can check that (16.6) is just a rewrite of (16.5): exercise 16.2 



p(xo) 



xo=f-'W 



|det/ f (^ )l 



(1 -dimensional) 
(<i-dimensional) . 



(16.9) 



For a deterministic, invertible flow x has only one preimage jco; allowing for mul- 
tiple preimages also takes account of noninvertible mappings such as the 'stretch 
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& fold' maps of the interval, to be discussed briefly in example 16.1, and in more 
detail in sect. 11.3. 



We shall refer to the integral operator with singular kernel (16.6) as the Perron- 
Frobenius operator. exercise 1 6.3 

example 23.7 



£ t (x,y) = d(x-f t (y)) . 



(16.10) 



The PeiTon-Frobenius operator assembles the density p(x, t) at time t by going 
back in time to the density p(xo, 0) at time t = 0. If you do not like the word "ker- 
nel" you might prefer to think of £!(x, y) as a matrix with indices x, y, and index 
summation in matrix multiplication replaced by an integral over y, {£? o p) (x) = 
J dyX,'(x, y)p(y) . In the next example Perron-Frobenius operator is a matrix, and remark 19.4 
(16.14) illustrates a matrix approximation to the Perron-Frobenius operator. 



Example 16.1 Perron-Frobenius operator for a piecewise-linear map: Consider 
the expanding I -dimensional map f(x) of figure 16.3, a piecewise-linear 2-branch map 
with slopes Ao > 1 and Ai = -Ao/(Ao - 1) < -1 : exercise 16.7 



fix) 



fo(x) - Aqx , 
/i(j0 = Ai(1-jc), 



xeMo = [0,1 /A ) 
x e Mi =(l/Ao,l]. 



(16.11) 



Both f(Mo) and f(Mi) map onto the entire unit interval M = [0, 1]. We shall refer to 
any unimodal map whose critical point maps onto the "left" unstable fixed point xo as 
the "Ulam" map. Assume a piecewise constant density 



p(x) = 



po ifx e Mo 
pi ifx e Mi 



(16.12) 



As can be easily checked using (16.9), the Perron-Frobenius operator acts on this 
piecewise constant function as a [2x2] Markov matrix L with matrix elements exercise 16.1 

exercise 16.5 



l 

|Aol 



I 

|Ai| 



|Aol |Ail 



(16.13) 



stretching both p and p\ over the whole unit interval A. In this example the density is 
constant after one iteration, so L has only a unit eigenvalue e s ° = 1/|A | + l/|Ai| = 1, 
with constant density eigenvector p = p\. The quantities 1/|A |, l/|Ai| are, respec- 
tively, the fractions of state space taken up by the \Mo\, \M\\ intervals. This simple 
explicit matrix representation of the Perron-Frobenius operator is a consequence of the 
piecewise linearity of f, and the restriction of the densities p to the space of piece- 
wise constant functions. The example gives a flavor of the enterprize upon which we 
are about to embark in this book, but the full story is much subtler: in general, there 
will exist no such finite-dimensional representation for the Perron-Frobenius operator, 
(continued in example 1 7.4) 



fast track: 

sect. 1 6.4, p. 334 
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16.3 Why not just leave it to a computer? 



Another subtlety in the [dynamical systems ] theory is that 
topological and measure-theoretic concepts of genericity 
lead to different results. 



— John Guckenheimer 



(R. Artuso and P. Cvitanovic) 



To a student with a practical bent the above Example 16.1 suggests a strategy for 
constructing evolution operators for smooth maps, as limits of partitions of state 
space into regions Mi, with a piecewise-linear approximations to the dynamics 
in each region, but that would be too naive; much of the physically interesting 
spectrum would be missed. As we shall see, the choice of function space for p is 
crucial, and the physically motivated choice is a space of smooth functions, rather 
than the space of piecewise constant functions. 

All of the insight gained in this chapter and in what is to follow is nothing but 
an elegant way of thinking of the evolution operator, as a matrix (this point of 
view will be further elaborated in chapter 23). There are many textbook methods 
of approximating an operator £, by sequences of finite matrix approximations 
but in what follows the great achievement will be that we shall avoid construct- 
ing any matrix approximation to £, altogether. Why a new method? Why not 
just run it on a computer, as many do with such relish in diagonalizing quantum 
Hamiltonians? 

The simplest possible way of introducing a state space discretization, fig- 
ure 16.4, is to partition the state space M with a non-overlapping collection of 
sets Mj, i - I,. . . ,N, and to consider densities (16.2) piecewise constant on each 



where Xi( x ) 1S the characteristic function (16.1) of the set Mi. This piecewise 
constant density is a coarse grained presentation of a fine grained density p{x), 
with (16.2) 



The Perron-Frobenius operator does not preserve the piecewise constant form, but 
we may reapply coarse graining to the evolved measure 



Mi. 





dxp(x). 




I dx(£o p)(x) 
Mi 




dyd(x-f(y)) , 
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B'RUTO IWSEHSlTlVO METHOD: 



Figure 16.4: State space discretization approach to 




computing averages. 

exact sp«dtwn J», ^<j^'i 's, °" 



or 



Pi = Zj PJ- 



7=1 



\Mj\ 



In this way 



Uj = 



\Minf- l (Mj) 
IMi 



P =pL 



(16.14) 



is a matrix approximation to the Perron-Frobenius operator, and its leading left 
eigenvector is a piecewise constant approximation to the invariant measure. remark 16.3 

The problem with such state space discretization approaches is that they are 
blind, the grid knows not what parts of the state space are more or less important. 
This observation motivated the development of the invariant partitions of chaotic 
systems undertaken in chapter 1 1 , we exploited the intrinsic topology of a flow to 
give us both an invariant partition of the state space and a measure of the partition 
volumes, in the spirit of figure 3.7. 

Furthermore, a piecewise constant p belongs to an unphysical function space, 
and with such approximations one is plagued by numerical artifacts such as spu- 
rious eigenvalues. In chapter 23 we shall employ a more refined approach to 
extracting spectra, by expanding the initial and final densities p, p' in some basis 
ipo, (p\, ip%, ■ ■ ■ (orthogonal polynomials, let us say), and replacing £.{y,x) by its 
ip a basis representation L aj g = {(p a \£\(fp). The art is then the subtle art of finding 
a "good" basis for which finite truncations of ~L a p give accurate estimates of the 
eigenvalues of chapter 23 



Regardless of how sophisticated the choice of basis might be, the basic prob- 
lem cannot be avoided - as illustrated by the natural measure for the Henon map 
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(3.17) sketched in figure 16.5, eigenfunctions of £, are complicated, singular func- 
tions concentrated on fractal sets, and in general cannot be represented by a nice 
basis set of smooth functions. We shall resort to matrix representations of X and 
the (f a basis approach only insofar this helps us prove that the spectrum that we 
compute is indeed the correct one, and that finite periodic orbit truncations do 
converge. 

in depth: 
chapter 1 , p. 3 

16.4 Invariant measures 

A stationary or invariant density is a density left unchanged by the flow 

p(x, t) = p(x, 0) = p(x) . (16. 15) 

Conversely, if such a density exists, the transformation f'(x) is said to be measure- 
preserving. As we are given deterministic dynamics and our goal is the compu- 
tation of asymptotic averages of observables, our task is to identify interesting 
invariant measures for a given f'{x). Invariant measures remain unaffected by dy- 
namics, so they are fixed points (in the infinite-dimensional function space of p 
densities) of the Perron-Frobenius operator (16.10), with the unit eigenvalue: exercise 16.3 

£<p(x)= f dyS(x-f t (y))p(y)=p(x). (16.16) 
JM 

We will construct explicitly such eigenfunction for the piecewise linear map in 
example 17.4, with p(y) = const and eigenvalue 1. In general, depending on the 
choice of f{x) and the function space for p(x), there may be no, one, or many 
solutions of the eigenfunction condition (16.16). For instance, a singular measure 
dp(x) = 5{x - x q )dx concentrated on an equilibrium point x q = f{x q ), or any 
linear combination of such measures, each concentrated on a different equilib- 
rium point, is stationary. There are thus infinitely many stationary measures that 
can be constructed. Almost all of them are unnatural in the sense that the slightest 
perturbation will destroy them. 

From a physical point of view, there is no way to prepare initial densities 
which are singular, so we shall focus on measures which are limits of transforma- 
tions experienced by an initial smooth distribution p(x) under the action of /, 

poW-lim f dy6(x-f(y))p(y,0), f dyp{y,0) = \. (16.17) 

Intuitively, the "natural" measure should be the measure that is the least sensitive 
to the (in practice unavoidable) external noise, no matter how weak, or round-off 
errors in a numerical computation. 
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16.4.1 Natural measure 



Huang: Chen-Ning, do you think ergodic theory gives us 
useful insight into the foundation of statistical mechanics? 
Yang: I don't think so. 

— Kerson Huang, C.N. Yang interview 

In computer experiments, as the Henon example of figure 16.5, the long time evo- 
lution of many "typical" initial conditions leads to the same asymptotic distribu- 
tion. Hence the natural measure (also called equilibrium measure, SRB measure, 
Sinai-Bowen-Ruelle measure, physical measure, invariant density, natural density, 
or even "natural invariant") is defined as the limit exercise 16.8 

exercise 16.9 



p xn (y) = 



lim^co i dr 6(y - f T (xo)) flows 

(16.18) 



lim^oo i HfcJ 6(y - f k (x )) maps , 

where xq is a generic initial point. Generated by the action of /, the natural 
measure satisfies the stationarity condition (16.16) and is thus invariant by con- 
struction. 

Staring at an average over infinitely many Dirac deltas is not a prospect we 
cherish. From a computational point of view, the natural measure is the visitation 
frequency defined by coarse-graining, integrating (16.18) over the Mj region 

Aft- = lim-, (16.19) 

f— > oo t 

where t\ is the accumulated time that a trajectory of total duration t spends in the 
Mi region, with the initial point xq picked from some smooth density p(x). 

Let a = a{x) be any observable. In the mathematical literature a(x) is a func- 
tion belonging to some function space, for instance the space of integrable func- 
tions L , that associates to each point in state space a number or a set of numbers. 
In physical applications the observable a(x) is necessarily a smooth function. The 
observable reports on some property of the dynamical system. Several examples 
will be given in sect. 17.1. 

The space average of the observable a with respect to a measure p is given by 
the (i-dimensional integral over the state space At: 

(a) P = ] r f dxp(x)a(x) 

\Pm\Jm 

\Pm\ - I dx p(x) = mass in M . (16.20) 
JM 

For now we assume that the state space M has a finite dimension and a finite 
volume. By its construction, (a) p is a function(al) of p. For p = po natural measure 
we shall drop the subscript in the definition of the space average; (a) p = (a). 
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Inserting the right-hand-side of (16.18) into (16.20), we see that the natural 
measure corresponds to a time average of the observable a along a trajectory of 
the initial point xq, 



*xo 



Uml r 

t^oo t Jo 



dra(f(x Q )). 



(16.21) 



Analysis of the above asymptotic time limit is the central problem of ergodic 
theory. The Birkhoff ergodic theorem asserts that if an invariant measure p ex- 
ists, the limit a(xo) for the time average (16.21) exists for (almost) all initial xq. 
Still, Birkhoff theorem says nothing about the dependence on xq of time averages 
a~^ (or, equivalently, that the construction of natural measures (16.18) leads to a 
"single" density, independent of xq). This leads to one of the possible definitions 
of an ergodic evolution: / is ergodic if for any integrable observable a in (16.21) 
the limit function is constant. If a flow enjoys such a property the time averages 
coincide (apart from a set of p measure 0) with space averages 



lim 

t— >oo 



- f draif 
t Jo 



(x Q )) = (a) 



(16.22) 



remark 16.1 
appendix A 



For future reference, we note a further property that is stronger than ergodicity: 
if the space average of a product of any two variables decorrelates with time, section 22.3 



lim (a(x)b(f(x))) = (a) (b) , 



(16.23) 



the dynamical system is said to be mixing. The terminology may be understood 
better once we consider as the pair of observables in (16.23) characteristic func- 
tions of two sets and S: then (16.23) may be written as 



lim 

t— >oo 



p(S) 



so that the set S spreads "uniformly" over the whole state space as t increases. 
Mixing is a fundamental notion in characterizing statistical behavior for dynam- 
ical systems: suppose we start with an arbitrary smooth nonequilibrium distribu- 
tion p(x)v(x): the after time t the average of an observable a is given by 

I dxp(x)v(f(x))a(x) 
JM 

and this tends to the equilibrium average (a) p if / is mixing. 



Example 16.2 The Henon attractor natural measure: A numerical calculation of 
the natural measure (16.19) for the Henon attractor (3.17) is given by the histogram 
in figure 16.5. The state space is partitioned into many equal-size areas Mi, and the 
coarse grained measure (16.19) is computed by a long-time iteration of the Henon map, 
and represented by the height of the column over area Mi. What we see is a typical 
invariant measure - a complicated, singular function concentrated on a fractal set. 
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Figure 16.5: Natural measure (16.19) for the Henon 
map (3.17) strange attractor at parameter values 
(a,b) = (1.4,0.3). See figure 3.6 for a sketch of the 
attractor without the natural measure binning. (Cour- 
tesy of J. -P. Eckmann) 




0.4 



-1.5 -0.4 



remark A. 2 



If an invariant measure is quite singular -for instance a Dirac 5 concentrated 
on a fixed point or a cycle- it is most likely of no physical import. No smooth 
initial density will converge to this measure if its neighborhood is repelling. In 
practice the average (16.18) is problematic and often hard to control, as generic 
dynamical systems are neither uniformly hyperbolic nor structurally stable: it is 
not known whether even the simplest model of a strange attractor, the Henon 
attractor of figure 16.5, is "strange," or merely a transient to a very long stable 
cycle. exercise 6.4 

16.4.2 Determinism vs. stochasticity 

While dynamics can lead to very singular p's, in any physical setting we cannot 
do better than to measure p averaged over some region At,; the coarse-graining is 
not an approximation but a physical necessity. One is free to think of a measure 
as a probability density, as long as one keeps in mind the distinction between de- 
terministic and stochastic flows. In deterministic evolution the evolution kernels 
are not probabilistic; the density of trajectories is transported detenninistically . 
What this distinction means will became apparent later: for deterministic flows chapter 1 9 
our trace and determinant formulas will be exact, while for quantum and stochas- 
tic flows they will only be the leading saddle point (stationary phase, steepest 
descent) approximations. 

Clearly, while deceptively easy to define, measures spell trouble. The good 
news is that if you hang on, you will never need to compute them, at least not 
in this book. How so? The evolution operators to which we next turn, and the 
trace and determinant formulas to which they will lead us, will assign the correct 
weights to desired averages without recourse to any explicit computation of the 
coarse-grained measure Ap,-. 



16.5 Density evolution for infinitesimal times 

Consider the evolution of a smooth density p(x) = p(x, 0) under an infinitesimal 
step 5t, by expanding the action of £ St to linear order in 8r. 
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= I dx 6(y - x - 6tv(x)) p(x) 
JM 

p(y - dTvjy)) _ p(y) - 5TVj(y)dip(y) 
|det (l + Sr^)\ " 1 + SrZU ^'OO 

n 

P (x,<5t) = p(x,0) - 6t—(v(x) P (x,0)) ■ (16.24) 

ox 

Here we have used the infinitesimal form of the flow (2.7), the Dirac delta Jaco- exercise 4.1 
bian (16.9), and the In det = tr In relation. By the Einstein summation conven- 
tion, repeated indices imply summation, vfy)di = Y^=\ v i(y)di- Moving p(y, 0) to 
the left hand side and dividing by 8t, we discover that the rate of the deformation 
of p under the infinitesimal action of the Perron-Frobenius operator is nothing but 
the continuity equation for the density: 

d t p + d ■ {pv) = . (16.25) 

The family of Perron-Frobenius operators operators {-£ f } feR+ forms a semigroup 
parameterized by time 

(a) £° = / 

(b) £'£'' = £ t+t ' t, f > (semigroup property) . 

From (16.24), time evolution by an infinitesimal step 5t forward in time is gener- 
ated by 

ftp(x) = + lim —(jfT-Apix) = -di( Vi (x)p(x)) . (16.26) 

<Sr->0+ OT v ' 

We shall refer to 



J{ = -d ■ v + V Vi(x)di (16.27) 



as the time-evolution generator. If the flow is finite-dimensional and invertible, 
J{ is a generator of a full-fledged group. The left hand side of (16.26) is the 
definition of time derivative, so the evolution equation for p(x) is 

|^-^ljp(x) = 0. (16.28) 

The finite time Perron-Frobenius operator (16.10) can be formally expressed 
by exponentiating the time evolution generator as 

£' = e m . (16.29) 
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The generator 3K is reminiscent of the generator of translations. Indeed, for a con- 
stant velocity field dynamical evolution is nothing but a translation by (time x velocity): 

exercise 16.10 

e~ ,v ^a{x) = a(x-tv). (16.30) 



16.6 Liouville operator 



A case of special interest is the Hamiltonian or symplectic flow defined 
by Hamilton's equations of motion (7.1). A reader versed in quantum mechan- 
ics will have observed by now that with replacement — > —iH, where H is 
the quantum Hamiltonian operator, (16.28) looks rather like the time dependent 
Schrodinger equation, so this is the right moment to figure out what all this means 
for Hamiltonian flows. 

The Hamilton's evolution equations (7.1) for any time-independent quantity 
Q - Q(q> P) are given by 



dQ = dQd^ + dQdn = dHdQ _ dQdH 
dt dqi dt dpi dt dpi dqi dpi dqi ' 



where (pt, qO span the full state space, which for Hamiltonian flows we shall refer 
to as the phase space. As equations with this structure arise frequently for sym- 
plectic flows, it is convenient to introduce a notation for them, the Poisson bracket 



remark 16.4 



, A m OA dB OA dB 

dpt dq t dqi dp t 

In terms of Poisson brackets the time-evolution equation (16.3 1) takes the compact 
form 

*§.={H,Q}. (16.33) 
dt 

The discussion of sect. 16.5 applies to any deterministic flow. The full phase 
space flow velocity is x = v = (q,p), where the dot signifies time derivative. section 28.1 

If the density itself is a material invariant, combining 
d,l + vdl = 0. 
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and (16.25) we conclude that <3,v, = and detJ'(xo) = 1- An example of such 
incompressible flow is the Hamiltonian flow. For incompressible flows the con- 
tinuity equation (16.25) becomes a statement of conservation of the phase space 
volume (see sect. 7.3), or the Liouville theorem 

dtp + Vjdip = 0. (16.34) 

The symplectic structure of Hamilton's equations (7.1) implies that the flow 
is incompressible, <9,v, = 0, so for Hamiltonian flows the equation for p reduces to 
the continuity equation for the phase-space density: 

d t p + di(pvi) = , i = \,2...,D. (16.35) 

Consider the evolution of the phase-space density p of an ensemble of nonin- 
teracting particles; the particles are conserved, so 

d id d d\ 

—p(q,p,t) = \ —+q i —+p i —\p{q,p,t) = 0. 

Inserting Hamilton's equations (7.1) we obtain the Liouville equation, a special 

case of (16.28): remark 16.4 

—p(q, p, t) = -ftp{q, p, t) = {H,p(q, p, t)} , (16.36) 
dt 

where { , } is the Poisson bracket (16.32). The generator of the flow (16.27) is in 
this case a generator of infinitesimal symplectic transformations, 

d d dH d dH d 
dqi dpi opidqi dq t dpi 

For example, for separable Hamiltonians of form H = p 2 /2m+ V(q), the equations 
of motion are 

Pi . dV(q) 

<?, = —, Pi = — o — ■ (16.38) 
m dqi 

and the action of the generator exercise 16.1 1 

Ji = --—+diV{q)—. (16.39) 
m dqi dpi 

Looking back at (16.30) we see that the first term generates a translation in the 
configuration space, f(q, p) — > f(q - dt q, p), and the second generates acceler- 
ation by force dV(q) in the momentum space. They do not commute, hence the 
time integration is not trivial. 
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The time-evolution generator (16.27) for the case of symplectic flows is called 
the Liouville operator. You might have encountered it in statistical mechanics, 
while discussing what ergodicity means for 10 23 hard balls. Here its action will 
be very tangible; we shall apply the Liouville operator to systems as small as 1 or 
2 hard balls and to our surprise learn that this suffices to already get a bit of a grip 
on foundations of the nonequilibrium statistical mechanics. 

Resume 

In physically realistic settings the initial state of a system can be specified only to 
a finite precision. If the dynamics is chaotic, it is not possible to calculate the long 
time trajectory of a given initial point. Depending on the desired precision, and 
given a deterministic law of evolution, the state of the system can then be tracked 
for a finite time only. 

The study of long-time dynamics thus requires trading in the evolution of a 
single state space point for the evolution of a measure, or the density of repre- 
sentative points in state space, acted upon by an evolution operator. Essentially 
this means trading in nonlinear dynamical equations on a finite dimensional space 
x - (jci , %2 • • • Xd) for a linear equation on an infinite dimensional vector space of 
density functions p(x). For finite times and for maps such densities are evolved by 
the Perron-Frobenius operator, 

p(x,t) = (£ t op)(x), 

and in a differential formulation they satisfy the continuity equation: 

d t p + d ■ (pv) = . 

The most physical of stationary measures is the natural measure, a measure robust 
under perturbations by weak noise. 

For long times the dynamics is described in terms of stationary measures, i.e., 
fixed points of the appropriate evolution operators. Reformulated this way, clas- 
sical dynamics takes on a distinctly quantum-mechanical flavor. If the Lyapunov 
time (1.1), the time after which the notion of an individual deterministic trajectory 
loses meaning, is much shorter than the observation time, the "sharp" observables 
are those dual to time, the eigenvalues of evolution operators. This is very much 
the same situation as in quantum mechanics; as atomic time scales are so short, 
what is measured is the energy, the quantum-mechanical observable dual to the 
time. Both in classical and quantum mechanics one has a choice of implementing 
dynamical evolution on densities ("Schrddinger picture," sect. 16.5) or on observ- 
ables ("Heisenberg picture," sect. 17.2 and chapter 18). 

In what follows we shall find the second formulation more convenient, but the 
alternative is worth keeping in mind when posing and solving invariant density 
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problems. However, as classical evolution operators are not unitary, their eigen- 
function can be quite singular and difficult to work with. In what follows we 
shall learn how to avoid dealing with these eigenstates altogether. As a matter of 
fact, what follows will be a labor of radical deconstruction; after having argued so 
strenuously here that only smooth measures are "natural," we shall merrily pro- 
ceed to erect the whole edifice of our theory on periodic orbits, i.e., objects that 
are 5-functions in state space. The trick is that each comes with an interval, its 
neighborhood - periodic points only serve to pin these intervals, just as millime- 
ter markings on a measuring rod are used to partition a continuum into intervals. 

Commentary 

Remark 16.1 Ergodic theory: An overview of ergodic theory is outside the scope of 
this book: the interested reader may find it useful to consult refs. [16.1, 16.3, 16.4, 16.5]. 
The existence of time average (16.21) is the basic result of ergodic theory, known as the 
BirkhofF theorem, see for example refs. [16.1, 16.25], or the statement of theorem 7.3.1 
in ref. [16.12]. The natural measure (16.19) of sect. 16.4.1 is often referred to as the SRB 
or Sinai-Ruelle-Bowen measure [1.29, 1.28, 1.32]. 

There is much literature on explicit form of natural measure for special classes of 
1-dimensional maps [1.19, 16.14, 16.15] - J. M. Aguirregabiria [16.16], for example, 
discusses several families of maps with known smooth measure, and behavior of measure 
under smooth conjugacies. As no such explicit formulas exist for higher dimensions and 
general dynamical systems, we do not discuss such measures here. 

Remark 16.2 Time evolution as a Lie group: Time evolution of sect. 16.5 is an 
example of a 1-parameter Lie group. Consult, for example, chapter 2. of ref. [16.13] for a 
clear and pedagogical introduction to Lie groups of transformations. For a discussion of 
the bounded semigroups of page 358 see, for example, Marsden and Hughes [16.6]. 

Remark 16.3 Discretization of the Perron-Frobenius operator operator It is an old 

idea of Ulam [16.18] that such an approximation for the Perron-Frobenius operator is a 
meaningful one. The piecewise-linear approximation of the Perron-Frobenius operator 
(16.14) has been shown to reproduce the spectrum for expanding maps, once finer and 
finer Markov partitions are used [16.19, 16.23, 16.20]. The subtle point of choosing a 
state space partitioning for a "generic case" is discussed in ref. [16.21, 23.22]. 



Remark 16.4 The sign convention of the Poisson bracket: The Poisson bracket 
is antisymmetric in its arguments and there is a freedom to define it with either sign 
convention. When such freedom exists, it is certain that both conventions are in use and 
this is no exception. In some texts [1.8, 16.7] you will see the right hand side of (16.32) 
defined as {B,A} so that (16.33) is ^ = {Q,H}. Other equally reputable texts [16.24] 
employ the convention used here. Landau and Lifshitz [16.8] denote a Poisson bracket by 
[A, B], notation that we reserve here for the quantum-mechanical commutator. As long as 
one is consistent, there should be no problem. 
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Remark 16.5 "Anon it lives"? "Anon it lives" refers to a statue of King Leontes's wife, 
Hermione, who died in a fit of grief after he unjustly accused her of infidelity. Twenty 
years later, the servant Paulina shows Leontes this statue of Hermione. When he repents, 
the statue comes to life. Or perhaps Hermione actually lived and Paulina has kept her 
hidden all these years. The text of the play seems deliberately ambiguous. It is probably 
a parable for the resurrection of Christ. (John F. Gibson) 



Exercises 



16.1. Integrating over Dirac delta functions. Check the 
delta function integrals in 

(a) 1 dimension (16.7), 

fdxS(h(x)) = J] -J-, (16.40) 

J {x:h(x>0} Wl 

(b) and in d dimensions (16.8), h : R d -> R d , 



f 

J n; 1 



dx 6(h(x)) 



Im 



z 



dx 5(h(x)) 
1 



det 



dh{ 



W- 41) 



{x:h(x)=0) I" 1 " 1 dx 

where Mj are arbitrarily small regions enclosing 
the zeros Xj (with xj not on the boundary dMj). 
For a refresher on Jacobian determinants, read, for 
example, Stone and Goldbart Sect. 12.2.2. 

(c) The delta function can be approximated by a se- 
quence of Gaussians 



dx5(x)f(x) = lim I dx 



■f' 



]2ncr 



--fix) . 



Use this approximation to see whether the formal 
expression 

dx 5{x 2 ) 



Consider 



makes sense. 



16.2. Derivatives of Dirac delta functions. 

Using integration by parts, determine the value of 



dxb{x)6 i2 \y) = Yj m 



{x:y(x)=Q} 



Vy" 



ly'l {(y'Y 



+b\3 



(y"Y 



(y'Y 

16.44) 



16.3. 



(y'f (y'Y 

These formulas are useful for computing effects of weak 
noise on deterministic dynamics [16.9]. 

generates a semigroup. Check that the Perron- 
Frobenius operator has the semigroup property, 

dz£ h (y, z) £ fl (z, x) = £ h+ " (y, x) , t u t 2 >0.(l> 



jr. 

Jm 



As the flows in which we tend to be interested are in- 
vertible, the £'s that we will use often do form a group, 
with t\,h e R. 

16.4. Escape rate of the tent map. 

(a) Calculate by numerical experimentation the log of 
the fraction of trajectories remaining trapped in 
the interval [0, 1] for the tent map 

/(x) = a(l-2|x-0.5|) 

for several values of a. 

(b) Determine analytically the a dependence of the es- 
cape rate y(a). 

(c) Compare your results for (a) and (b). 

16.5. Invariant measure. We will compute the invariant 
measure for two different piecewise linear maps. 



dx 8'(y) 
dx^ 2) (y) 



where y = /(x) - jc (16.42) 
1 



Z 



[x:y(x)=0) 



,(y") 2 



[y'l I (y'Y 
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(a) Verify the matrix X representation (16. 13). 

(b) The maximum value of the first map is 1 . Com- 
pute an invariant measure for this map. 

(c) Compute the leading eigenvalue of X. for this map. 

(d) For this map there is an infinite number of in- 
variant measures, but only one of them will be 
found when one carries out a numerical simula- 
tion. Determine that measure, and explain why 
your choice is the natural measure for this map. 

(e) In the second map the maximum occurs at a = 
(3 - V5)/2 and the slopes are +( V5 + l)/2. Find 
the natural measure for this map. Show that it is 
piecewise linear and that the ratio of its two values 
is(V5 + l)/2. 

(medium difficulty) 

16.6. Escape rate for a flow conserving map. Adjust A , 
A! in (16.1 1) so that the gap between the intervals Aio, 
Mi vanishes. Show that the escape rate equals zero in 
this situation. 

16.7. Eigenvalues of the Perron-Frobenius operator for the 
skew full tent map. Show that for the skew full tent 
map 




the eigenvalues are available analytically, compute the 
first few. 

16.8. "Kissing disks"* (continuation of exercises 8.1 and 
8.2) Close off the escape by setting R — 2, and look 
in real time at the density of the Poincare section iter- 
ates for a trajectory with a randomly chosen initial con- 
dition. Does it look uniform? Should it be uniform? 
(Hint - phase-space volumes are preserved for Hamil- 
tonian flows by the Liouville theorem). Do you notice 
the trajectories that loiter near special regions of phase 
space for long times? These exemplify "intermittency," 
a bit of unpleasantness to which we shall return in chap- 
ter 24. 



16.9. Invariant measure for the Gauss map. 

the Gauss map: 



/(*) = 







[i] 



x±0 
x = 



Consider 



(16.47) 



where [ ] denotes the integer part. 

(a) Verify that the density 
1 1 



p(x) = 



log 2 1+ i 
is an invariant measure for the map. 
(b) Is it the natural measure? 

16.10. M as a generator of translations. Verify that for 
a constant velocity field the evolution generator in 
(16.30) is the generator of translations, 



/(*) = 



fo(x) = A x , 



x), 



16.11. 



*eMo = [0,l/Ao) 
xeMi =(1/A ,1]. 



e tVd -a{x) = a(x + tv) . 

Incompressible flows. Show that (16.9) implies that 
po(x) = 1 is an eigenfunction of a volume-preserving 
flow with eigenvalue sq = 0. In particular, this im- 
plies that the natural measure of hyperbolic and mixing 
Hamiltonian flows is uniform. Compare this results with 
^^'rramerical experiment of exercise 16.8. 
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